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Friction characteristics with respect to surface topographic orientation were investigated using surfaces of 
different materials and fabricated with grooves of different scales. Scratching friction tests were conducted 
using a nano-indentation-scratching system with the tip motion parallel or perpendicular to the groove 
orientation. Similar friction anisotropy trends were observed for all the surfaces studied, which are (1) under 
a light load and for surfaces with narrow grooves, the tip motion parallel to the grooves offers higher friction 
coefficients than does that perpendicular to them, (2) otherwise, equal or lower friction coefficients are 
found under this motion. The influences of groove size relative to the diameter of the mating tip (as a 
representative asperity), surface contact stiffness, contact area, and the characteristic stiction length are 
discussed. The appearance of this friction anisotropy is independent of material; however, the boundary and 
the point of trend transition depend on material properties. 
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Anisotropic surfaces are reported on the skins of wild animals, sandfish^ and snake^'^ for example, and are 
results of many fabrication and manufacturing processes, such as turning, grinding, honing, and shaping, 
etc. Friction anisotropy probed in different motion direction is widely observed in interactions of various 
surfaces, such as single crystal surfaces'*"^, quasicrystal surfaces^'^, organic molecular monolayers'"^ \ monolayer 
graphene^^, and textured surfaces^^'^'*. It even exists in the presence of lubrication^^"^^. Differences in surface 
potential corrugations and bulk plasticity with respect to lattice structures''^' distinct excitations of electron and 
phonon in quasicrystals^'^'^'"^\ surface stiffness differences caused by molecular chain alignments ^°'^^ or graphene 
puckering^^, and micro-hydrodynamic lifting and fluid side flows^^'^^'^^ have been analyzed in the course of 
understanding the nature of friction anisotropy. 

Friction anisotropy may involve complicated physical phenomena, but in many cases, surface topography 
plays a crucial rule and results in friction anisotropy with respect to topographic orientation. This type of friction 
anisotropy is gaining increasingly strong attention due to its importance to innovative designs of frictional 
interfaces; however, no systematic research has been committed so far to reveal the scientific nature of related 
friction variation. 

A grooved texture surface mating with a tip of a known size, as illustrated in Fig. la, leads to a contact interface 
formed by single or multiple surface grooves with a relatively smooth spherical asperity (the tip), depending on 
the scale of the surface textures (Figs, lb and Ic). The relative motion of the tip parallel and perpendicular to the 
grooves should reveal the sensitivity of friction to surface pattern orientation. The current study is committed 
with this setup, as a systematic comparative study, to investigate the friction anisotropy phenomena with respect 
to surface topographic characteristics and to understand the roles of pattern size, applied load, and materials. This 
study involves four materials, three scales of interactions, and a range of loads. Several key questions are addressed 
for the first time, such as whether friction anisotropy is a material related issue, would it change with topographic 
parameters, and is there a critical transition point of the friction anisotropy at which the effects of topographic 
orientations switch? The answers to these questions should shed light on friction control by means of surface 
patterning and texturing, similar to modifying the "landscape" of surface interaction to modulate friction at the 
atomic scale^l 

Results 

Friction coefficients of all samples were obtained with the same methods. Figure 2 shows the results for G 139 Si 
surfaces (bank width 139 nm, pitch 278 nm) and G4 WC surfaces (bank width 110 |im, pitch 220 |im). In the 
perpendicular direction, both the friction coefficient history and the real-time tip trajectory fluctuate simulta- 
neously following the texture profile. Although the period of the curves for the nano-scale pattern of G 139 Si 
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Figure 1 | Surface orientations and groove patterns, (a) Schematic of two characteristic surface orientations with respect to the relative motion 
directions (the drawing is not to scale). The bank width a and groove width h are shown, (b) SEM images of the G139 and G278 grooves in silicon surfaces 
(courtesy of LightSmyth Technologies, Inc.). (c) SEM images of the G1-G4 textures in 52100 steel surfaces, (d) Optical measurement results of the Gl 
52100 steel surface with topography (left) and cross section profile (right). The pitch, defined as a+ b, of the Gl texture is confirmed to be about 4 |am. 



samples are not distinguishable considering the small tip motion 
(4.0E-03 nm as Usted in Table SI of the Supplementary material) 
in the depth direction and large friction variation, that for the G4 WC 
samples is clearly the same as the groove pitch (Fig. 2a) of about 
220 |im, and the shapes of both curves in the parallel direction 
(Fig. 2b) are nearly flat. The tip motion trajectories in the two char- 
acteristic directions of the G139 Si surfaces are shown in nm to 
emphasize the details. In contrast to the G4 WC surfaces, the tip 
fluctuation on the G139 Si surface in the parallel direction is larger 
than that in the perpendicular direction. When the WC tip (root 
mean square (RMS) roughness =66.0 ±5.6 nm) travels in the parallel 
direction, stick- slip happened on the nearly atomic smooth continu- 
ous Si banks (RMS roughness =2.3 ±0.2 nm). The tip jump due to 
the release from stiction leads to the peaks in the tip trajectory. The 
discontinuous patterns in the perpendicular direction prohibit the 
tip from sticking to the grooved surface. However, no obvious stick 
and slip were observed on the G4 WC surface in the parallel direction 
due to the higher surface roughness (RMS roughness values are 
151.0±2.0 nm for the 304 stainless steel tip and 10.3±2.1 nm for 
the non-grooved WC surface). It is worth mentioning that the low- 
speed friction responses of the G4 surface in the perpendicular dir- 
ection is similar to a stick- slip behavior; however, the energy barrier 
imposed by the texture ridges, instead of stiction, is the source of 
friction variation. 

The friction coefficients were obtained by averaging data over a 
220 |im traverse distance for at least 3 repeats. Figure 2c shows that 
the average friction coefficients are direction dependent, with higher 
values in the parallel direction for samples with narrower grooves, 
such as G139, and equal or higher values in the perpendicular 



direction for samples of wider grooves of G4. The large friction error 
bar in the parallel direction of the G4 samples is due to the location 
that the tip contact was made, e.g. on the bank or at the corner edges 
of a groove, and higher friction is associated with the contact on the 
bank as if there were no textures. 

Friction coefficients of all samples are given in Fig. 3. The average 
friction coefficients for three types of surfaces, i.e., copper alloy (Cu), 
52100 steel and tungsten carbide (WC), with Gl to G4 textures in two 
directions are plotted against the relative pattern size, b/R, groove 
width b mating over tip radius R, in Fig. 3a. This figure reveals that 
the comparison of the average friction coefficients from the parallel 
and perpendicular directions of motion depends greatly on the 
groove size with respect to the tip size. The average friction coeffi- 
cients for the Gl surfaces plotted at the left end of Fig. 3a show higher 
values in the parallel direction; while those for the G4 surfaces, at the 
right end, are higher in the perpendicular direction. Variation of the 
average friction coefficient points out a relative size, b/R=(b/R)c, for 
the transition of the friction trend; assuming a perfect rectangle 
profile without deformation, the friction coefficient and topographic 
parameters, such as groove width b and tip radius R, are related as 

Mave~M|l + — l] where fi^ve is the average friction 

coefficient, and fi is the friction coefficient related to a flat surface. It 
obviously indicates that the average friction coefficient increases with 
b/R, which explains the geometric meaning of the transition. This 
observation is in agreement with the higher friction reported at the 
edges of atomic-scale surface steps due to deeper potential bar- 
riers^^'^^ Moreover, the value of (b/R)^ changes with material effective 
modulus, suggesting that surface deformation, or issues related to 
surface stiffness, might have played a role in the transition. 
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Figure 2 | Friction coefficient histories and real-time tip trajectories, (a) Friction coefficient history curves and real-time tip trajectories in the 
perpendicular direction. The groove surface profile of G4 WC is plotted for reference, (b) Friction coefficient history curves and real-time tip trajectories 
in the parallel direction. Note that the tip trajectory fluctuations are around 5 to 10 nano meters in the G 139 Si surfaces, much smaller than the groove size, 
(c) Average friction coefficients in the parallel (//) and perpendicular ( _L ) directions of a G139 Si and a G4 WC surface. The experiment details are 
following: applied load: 5 mN, tip speed: 1 |am/s. The G139 Si sample (bank width 139 nm, pitch 278 nm) was mated with a WC ball (R=803± 12 jam, 
RMS roughness =66 ±5.6 nm); the G4WC sample (band width 110 |am, pitch of 220 |am) was mated with a 304 stainless steel ball (R= 800 ±3 jam, RMS 
roughness=151±2 nm). 



Figure 3b plots the average friction coefficients with respect to 
normalized applied load defined as load over the critical load at yield, 
or a failure threshold of the corresponding smooth surface. The 
friction anisotropy phenomenon is revealed for each material pair 
under various loads. Note that the difference of friction coefficients 
reduces and the anisotropy vanishes, or the trend even reverses at a 
certain point, as load increases, which suggests boundaries of such 
friction anisotropy. The dependence of the position of an anisotropy 
boundary on load implies the effect of surface deformation, and 
hence the influence of modulus and hardness. The smoother tip 
trajectory on the G2 Cu sample under the 100 mN load reveals that 
the texture was smoothened (Fig. 3c) as a result of significant 
deformation. The textures in the 52100 steel surfaces were still in a 
good shape even after experiencing the maximum applied load of the 
test range, and for this reason the friction anisotropy remains. A 
higher load is needed to explicitly capture the boundaries of the 
friction anisotropy of this material. The decreasing amplitudes of 
the tip trajectories of the 52100 steel surfaces explain the reduction 
of the friction anisotropy. The average friction coefficients associated 
with the G2 WC surfaces in two directions merge at the normalized 
load of about 0.5; then they deviate again but in the opposite dir- 
ection. This is attributed to the deformation occurred in the relatively 
more compliant stainless steel tip. Rubbing while deforming should 
increase the energy barrier for the motion along the perpendicular 
direction, as shown by the sharpened tip trajectories. 

Macroscopic friction tests was conducted to confirm the above 
mentioned findings using a universal material tester (CETR, UMT- 
2) with a 6.35 mm diameter Cr ball indenting and sliding on the G2 
WC surfaces under 196 mN load with sliding speed of 50 jim/s. The 
friction coefficient curve in the perpendicular direction has a period 
of about 16 jim, consistent with the groove size. The average friction 
coefficients in the parallel and perpendicular directions were 
0.0422±0.0007 and 0.0331 ±0.0022, respectively. The result are in 



agreement with the data at the left-hand sides of the WC plots in 
Fig. 3a and 3b before the transition. Thus, friction anisotropy with 
respect to topographic orientation is confirmed by these macroscopic 
friction tests, with higher friction in the parallel direction for hIR < 
(b/R)c and a normalized load (0.16) lower than the critical load ratio 
for the friction transition to occur (0.5). 

Discussion 

The friction enhancement in the perpendicular directions for sur- 
faces with wider groove patterns as a result of higher energy barriers 
is believed to lead to the transition of friction anisotropy. The reason 
why the friction coefficients in the parallel direction cases are higher 
than those in the perpendicular direction cases for surfaces with 
narrower patterns deserves an attention. In this work, moisture lub- 
rication is not likely due to low humidity. The contact materials are 
dissimilar materials, so that commensurability and mismatch of the 
crystalline directions are not applicable. Plastic deformation and 
plowing are not involved under light loads. The surfaces were well 
cleaned and the tracks were not re-used, consequently no obvious 
third body or wear debris problems should be considered. 

The strong orientation of grooves should lead to anisotropic sur- 
face stiffness. Finite element simulations were conducted using 
SIMULIA Abaqus FEA, and the results show that the surface contact 
stiffness, inversely proportional to surface deformation, is larger 
when the tip motion is parallel to the grooves than that when the 
tip motion is perpendicular to the grooves. Among the microscale 
G1-G4 52100 steel samples, for example, the maximum ratio of the 
surface stiffness (// over ±) is about 1.13 for the Gl textured surfaces, 
while the ratios for other 52100 samples decrease with increased 
bank widths, which should finally reach unit for very wide grooves. 
The surface stiffness ratios are normalized by their maximum value 
among all textured surfaces and shown as circles in Fig. 4. The same 
trend is observed for other samples although that for the nanoscale 
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Figure 3 | Friction characteristics from size-dependent and load-dependent tests, (a) Friction anisotropy under the constant load of 5 mN with respect 
to the relative pattern size, h/R, corresponding to the groove bank sizes of 2, 8, 25 and 110 jam. (b) Friction anisotropy of the same grooves (8 |am) with 
respect to normalized load, defined as applied load over the critical load of yielding or failure, corresponding to the applied loads of 5, 10, 30 and 100 mN. 
(c) Real-time tip trajectories from the load-dependent tests. Note that the curves are separated with certain distance for clarity. The mating tips are 304 
stainless steel balls (R=800±3 |am, RMS roughness=151±2 nm), and the scan speed was 1 |am/s for all of the above experiments. 



textured Si surfaces follows a different line. Trend discontinuity is 
observed in between the Gl and G278 samples as a result of pattern 
sizes in two different scales, which causes distinct number of banks in 
contact for these surfaces. When several texture banks are in contact, 
the load supported by each of them is quite different; therefore, the 
stiffness of each texture is different, and the overall surface stiffness is 
the combined action of all individual texture. However, the texture 
bank right underneath the tip takes the majority of the load. As a 
result, this bank has the maximum surface contact stiffness ratio (see 
Supplementary Fig. S2). Even with the surface stiffness impact exam- 
ined, the direct correlation between surface stiffness and friction for 
the material group in this study is not as clear as that for compliant 
materials, graphene for example. 

The evolution of contact areas follows the same manner as friction 
variations, i.e., continuously smooth in the parallel direction but 
fluctuating in the perpendicular direction. The degree of fluctuation 
depends on relative groove size h/R; slight change is associated with 
the case of a small h/R, while a larger fluctuation with the case of a 
larger h/R (see Fig. S3 in Supplementary Information). Since friction 
anisotropy holds even on surfaces with h/R as small as 1.7X10"^ 
(G139), the effect of this contact area difference in the two directions 
may be trivial. 



The critical stiction lengths in the two orthogonal directions 
obtained from the micro -slip analysis at the local contacts are also 
believed to be responsible for friction anisotropy^^. Compared with 
the discontinuous shorter contacts in the perpendicular motion 
cases, a longer stiction length is associated with the parallel cases 
prior to macro -slip initiation. At low- speed scratching, this means 
continuous adhesion in the contact junction when motion starts, and 
thus more energy has to be dissipated^^ and a high friction is expected 
in the cases of low topography- induced barrier. From the peeling 
of a flexible plate bonded to a rigid substrate with an incision - 
patterned thin adhesive, Chung^^ reported that longitudinal incisions 
demanded higher loads to propagate the cracks than lateral incisions 
did, which demonstrates the same trend as the friction anisotropy of 
surface with narrower grooves. Different crack propagation speeds 
and fracture contours from elliptical contacts in the shear mode 
along two directions were reported by Das et al.^^ which are analog- 
ous to friction anisotropy with respect to topography. 

A static friction model developed by Chen and Wang^^ for solving 
the contact between a sphere and a rough surface was employed to 
analyze the maximum stiction length, which is the longest continu- 
ous length in the stick area measured along the direction of the ten- 
dency of motion (see Fig. S4 in the Supplementary Information). If 
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Figure 4 | Relative friction anisotropy (COF change over COF _L) 
related to surface stiffness ratio (// over ±) and stiction length ratio 
(// over ±) in regard to surface orientations and energy-barrier induced 
COF change in the perpendicular direction. The results of each effect are 
normalized by its maximum value. The domain is classified as two groups: 
surfaces with nanoscale textures, G139 and G278, and other surfaces with 
microscale textures, G1-G4. In each region, the possible mechanisms of 
friction anisotropy are listed in the order of impact. 

the shear traction is less than the local shear strength, the local 
contact area is in stick, otherwise, in slip. The ratios of the maximum 
stiction length (// over ±) are normalized by their maximum value 
among all textured surfaces and shown as squares in Fig. 4. Among 
the microscale textured 52100 steel samples (G1-G4), the largest 
ratio before normalization is 3.41 for the Gl textured surface, while 
that of the rests are all unit because the contact widths are less than 
the bank widths for these cases. On the other hand, the stiction length 
of the nanoscale textured surfaces, G139, is dramatically reduced by 
the discontinuous patterns in the perpendicular orientation. 

The variations of friction, surface stiffness, and stiction length with 
respect to surface orientations and the COF enhancement due to 
energy barrier in the perpendicular direction were shown altogether 
in the Fig. 4. The friction anisotropy is presented as the relative 
change of COF over the COF in the perpendicular direction; the 
stiction length and surface stiffness results are also displayed as ratios 
(// over ±); the energy barrier induced COF enhancement is esti- 
mated by — ^~[4(^/i^)~^-l]~^ where \i and \iave are the 

friction coefficient associated with a smooth surface and that with 
a textured surface in perpendicular direction, respectively, and h/R is 
the relative grove size. The impacts of the surface stiffness, stiction, 
and engery barrier are normalized by their maximum values in Fig. 4. 
The result domain can be clearly divided into two regions, one for the 
nanoscale grooves (left) and the other for the microscale textures 
(right). It seems that the stiction length difference dominates the 
COF variations in the former for G139 and G278, while the effect 
of surface stiffness plays a secondary role and the groove-induced 
energy barrier in perpendicular direction is negligible because their 
patterns are the smallest. In the region of microscale textured sur- 
faces, G1-G4, the energy barrier effect is the strongest. The larger the 
relative pattern, the steeper the barrier; and this observation is sup- 
ported by the friction trend reverse for the G4 samples. The effects of 
stiction length and surface stiffness are both active at this range, but 
they become less important as the surface patterns become larger. 
The effect of surface stiffness seems have larger influence zone com- 
pared with that of stiction length because the surface stiffness ratios 
are still distinguishable among G2-G4 samples. 



To summarize, the exploration of this work covers a wide range of 
materials; similar friction anisotropy was observed for all the surfaces 
of the materials studied, and the same trend of transition of friction 
anisotropy with respect to topography was observed. The analyses 
reveal that the surface contact stiffness, the energy barrier in the 
perpendicular direction and the stiction length in the direction of 
motion have notable impacts on the friction anisotropy phenom- 
enon, and the contribution from each is closely related to the relative 
pattern size. 

Methods 

Groove fabrication and characterization. A series of groove- shape textures with 
bank width a of about 2 |im, 8 |am, 25 |am, and 110 |am, named as Gl to G4 
respectively, were fabricated on polished copper alloys (Cu), 52100 steel, and tungsten 
carbide (WC) surfaces by laser ablation with a pico-second laser system (RAPID 
Lumera Laser). Silicon stamps (from LightSmyth Technologies, Inc.) with nano 
grooves of 139 nm and 278 nm, labeled as G 139 and G278, were also available for 
study. Each sample is named based on its material and groove size, such as Gl WC for 
tungsten carbide surfaces with 2 |im grooves, G 139 Si for silicon surfaces with 
139 nm grooves. Figure lb and Ic show the sample groove designs and 
measurements. 

The grooved surfaces were examined by using a white-light interferometer (zygo, 
NewView 7300) and a scanning electron microscopy (SEM) (Hitachi, S-4500) prior to 
the friction tests. The coverage ratio, defined as bank width a over pitch a+ ^, where h 
is the groove size, is about 50% for all samples. The aspect ratio, which is defined as 
bank width a over groove height /z, ranges from 1 to 22 (shallow features). Geometric 
parameters of WC samples can be found in Supplementary Table SI online. The 
geometries of the other two materials are largely the same with very small variations 
due to the fact that laser ablation is material-dependent. The laser ablation method is 
capable of creating complicated patterns on various materials; however, the minimal 
feature size is constrained by the lens in use and tiny paths are generated on the 
bombarded surface, as shown in the SEM image of the Gl surface. The Gl profile 
measured with a white-light interferometer shows a sinusoidal-like cross section 
instead of the targeted perfectly rectangles. Therefore its real bank width at the top is 
less than the targeted value of 2 |am. However, the pitches were precisely controlled to 
be 4 |j,m for the Gl texture, and so were the others to be 16 |im, 50 |im and 220 |j,m 
for G2-G4 textures, respectively. The amplitude of tip trajectory A shows the tip 
motion in the height/depth directions and can be calculated from the geometric 
relationship between a spherical tip and a groove, assuming perfect rectangle profile 
without deformation. The depths of all the grooves are deep enough to prevent the tip 
from contacting the groove bottom land. 

Mating tips. The mating tips were made by rigidly gluing WC balls (to mate with Si 
surface) and 304 stainless steel balls (for the rests) onto precisely machined tip holders 
fitting to the test equipment to be described later. The white-light interferometry 
measurements reveal that the true radius and the surface root mean square (RMS) 
roughness are 800 ±3 |am and 151 ±2 nm for the 304 stainless steel balls, and 
803 ± 12 |im and 66 ±5.6 nm for the WC balls, respectively. 

Friction measurements. A nano-indentation-scratching system (NanoTest, Micro 
Materials Ltd, UK) was used to measure friction by slowly sliding a spherical indentor 
tip on a grooved sample surface in a controlled environment (about 26°C and 30% 
humidity). Each specimen was carefully aligned with the grooves either parallel or 
perpendicular to the indentor motion direction, as illustrated in Fig. la, with the 
assistance of the attached microscope. In the tests, normal loads of 5, 10, 30 and 
100 mN were used to investigate the load effect, and 5 mN was used for all other tests. 
The scanning speed was 1 |im/s. The indenter penetration depth, scan distance, and 
the normal and lateral forces were recorded during each test after reaching the steady 
state. Friction coefficient (or COF) is obtained as the ratio of the recorded lateral force 
over the applied normal load. At least three measurements (scratches) were taken at 
different locations under each operating condition, and the results were averaged. For 
the surfaces with large textures, such as G4, the tests in the parallel direction were 
conducted on the banks and over the grooves for 3 times each. A pin-on-disk tester 
was also used to measure macroscopic friction at the same environment. 
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